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Abstract: Broadband homodyne detection of the light transmitted by 
a Fabry-Perot cavity containing a strongly-coupled ^^^Cs atom is used 
to probe the dynamic optical response in a regime where semiclassical 
theory predicts bistability but strong quantum corrections should apply. 
While quantum fluctuations destabilize true equilibrium bistability, our 
observations confirm the existence of metastable states with finite lifetimes 
and a hysteretic response is apparent when the optical drive is modulated on 
comparable timescales. Our experiment elucidates remnant semiclassical 
behavior in the attojoule (~ 10 photon) regime of single-atom cavity QED, 
of potential significance for ultra-low power photonic signal processing. 
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Research over the past decade in single-atom cavity quantum electrodynamics (cQED) has 
largely focused on the generation of non-classical states of light and on the development of 
devices for quantum information protocols IB E El • However single-atom cQED with strong 
coupling is also a natural context in which to study the interplay of nonlinear mean-field dynam- 
ics and quantum fluctuations |i4|[5j|6l[7i|, which is a topic of equally broad fundamental interest 
and with potentially greater near-term relevance for information technology |8|. Dynamical 
systems theory |9| provides key methods and insights for nonequilibrium statistical mechanics, 
as well as essential tools for analyzing and engineering nonlinear phenomena. Substantial work 
will be required to generalize these methods to encompass quantum coherence and fluctuations 
for applications in the emerging disciplines of nanoscale and quantum engineering. Absent a 
systematic understanding of quantum stochastic nonlinear phenomena in the strong coupling 
regime, intuitive connections to semiclassical theory provide important guidance for analyzing 
and engineering the behavior of open quantum systems Lfi, J0|. Here we show that the optical 
response of a strongly driven single-atom cQED system clearly displays competing influences 
of semiclassical bistability and quantum fluctuations, in quantitative agreement with theory. 

Our experiment utilizes a single gas-phase ^^^Cs atom as the nonlinear medium in a Fabry- 
Perot optical resonator, but we note that the same physics should also be relevant to quan- 
tum nonlinear dynamics in strongly coupled nanophotonic systems ITT] ^| . We study the dy- 
namic input-output properties of the atom-cavity system in a parameter regime that semiclassi- 




Fig. 1. (a) Steady state intracavity mode amplitude a from the semiclassical Maxwell- 
Bloch Equations (MBEs) as a function of drive amplitude E for the experimental cQED 
system with {0, A} = {— 1.1, .Tj/c mode-drive and atom-drive detunings. Blue (green) in- 
dicates dynamically stable (unstable) solutions, predicting true amplitude bistability for 
drives in the interval E = [1.5, 2.3] K". (b) Wigner quasi-probability functions of the steady 
state cavity field from the quantum model as a function of drive amplitude for the same 
parameters as (a), (c) Sample trace of amplitude quadrature homodyne measurement of 
the transmitted field during an atom transit with the drive turned on at ^ = 2^s and held 
Sit E = 2.6k. The slight decrease in signal variance between 2 and 14/is is likely due to 
a gradual decrease in the coupling rate g from a initial, near-maximal value as the atom 
moves through the cavity mode. 

cally (13] would be expected to exhibit absorptive bistability with a hysteresis loop fT4\ suitable 
for attojoule (~ 10 photon) optical switching |15|. Theoretical studies have shown B1[T6l[6ll 
however that quantum fluctuations of the cavity field and atomic dipole should induce sponta- 
neous switching between low- and high-transmission states, destroying true optical bistability, 
when the energy separation of the semiclassical dynamical attractors reaches the few-photon 
scale. Previous experimental studies have proven the need to utilize fully quantum models for 
predicting the steady-state optical response in single-atom cQED with strong driving 1 17 |. Here 
we take an important step further by recording dynamic signals that reaffirm the quantum me- 
chanical model, but nonetheless reveal remnant behavior reminiscent of the semiclassical pic- 
ture even in the deep quantum regime. As research on nanophotonic logic devices pushes to- 
wards attojoule- scale switching energies (TSl [191 0, functionalizing such remnant behaviors 
will become a definitive challenge. 

The Max well-B loch Equations (MBEs) represent a semiclassical mean-field approximation 
to the cQED master equation (see Appendix [BJ in the limit of weak coupling and large atom 
number |[T3]|20l and have been used in research on optical bistability with atomic ensembles as 
the intracavity nonlinear medium fJT, 22, 6|. In single-atom cavity QED with strong coupling 
and saturating driving fields, the MBEs should not apply a priori as they ignore correlations 
between the atom and intra-cavity photons. However, the MBEs do retain some relevance as a 
projection of the quantum master equation onto the sub-manifold of semiclassical atom-field 
states |23 1 and numerical studies have shown 1 6 1 that solutions of the full quantum model often 
mimic those of the MBEs qualitatively. 

The MBEs predict absorptive optical bistability for the parameters of our experiment, as 
shown in Fig.[TJa). The vertical axis represents input drive amplitude (with ^/TJpW per unit on 
the vertical axis), and the solution curve passes through horizontal planar coordinates propo- 
tional to the real and imaginary parts of the output steady-state complex field amplitude(s). The 
projection of the curve onto the bottom face of the coordinate box indicates the phases of the 
steady-state solutions (relative to the optical phase of the laser drive), while the projection onto 



the back face illustrates the input-output relation for amplitude-quadrature homodyne detec- 
tion. In this representation, a characteristic 'S-curve' appears with two stable solutions and one 
unstable solution co-existing for E in the range of [1.5, 2.3] /c. Using identical parameters, the 
steady- state solution of the quantum master equation can be determined for each value of the 
drive strength. Contours of the corresponding intracavity field Wigner functions (tracing over 
the atomic states) are displayed for several drive strengths in Fig.[TJb). A double-peaked struc- 
ture emerges along the amplitude quadrature in a range of E/k similar to the bistable region 
of the MBEs. The bimodal steady states correspond to an incoherent mixture of two states: 
a weakly excited atom and low-amplitude field state, and a fully saturated atom and a high- 
amplitude field state. Thus, the two peaks of the distribution may be qualitatively associated 
with the low- and high- amplitude branches of the semiclassical absorptive bistability curve, but 
in the quantum model neither is truly stable | 4 |. It can be seen that in any single trial the cavity 
field (and thus the output power) spontaneously switches between low and high states, as in 
the experimental data of Fig. [TJc) (described below). This stochastic switching is a dramatic 
consequence of quantum fluctuations in this non-linear, attojoule- scale optical system. 

In our experiment ||7l[24l, laser cooled ^^^Cs atoms are dropped into a Fabry-Perot optical 
resonator (length 27/im, 10cm radius of curvature mirrors, field decay rate k/2tc = 9.3MHz) 
supporting a circularly-polarized, 852nm TEMqo mode actively frequency- stabilized relative 
to the {6Sii2-,E = A^mp = +4) (6/3/2,^ = S^mp = +5) atomic cycling transition (with 
dipole decay rate Y±/27l = 2.6MHz). As an atom falls through the cavity mode it experiences 
a position-dependent coupling rate g (maximum go/^n =56.8MHz at the cavity anti-nodes) 
that can be monitored via the transmission of a weak and detuned circularly -polarized optical 
probe. Once a strongly coupled atom is detected, the probe power and detuning are adjusted for 
optical homodyne detection and 200MS/s data acquisition (first arrow at 2/is on the time axis 
in Fig.[TJc). Fig.[TJc) depicts a representative signal (all data has been post-filtered at 20MHz 
bandwidth for clarity); with the drive amplitude held at a fixed value of £^ = 2.6^*, the amplitude 
quadrature of the transmitted field fluctuates with a large variance until 14/is (red arrow) when 
the atom is abruptly lost and and the measured transmission settles to the shot noise- variance 
signal with intermediate mean amplitude expected for our cavity when empty (perturbed only 
by slight ~l-10kHz mechanical instabilities). See Appendix [A| for more explanation of the 
apparatus. 

Zooming in on segments of the trace in Fig. [TJc) reveals significant dynamics on /is time 
scales. While the signals after the atom has left the cavity have white noise, Gaussian statistics 
as expected (e.g. the 14 — 19/is segment from Fig.[TJc) is shown in Fig.[2ja)) single shot signals 
observed when the drive is first turned on are more volatile (the 2 — 7/is segment from Fig.[TJc) 
is shown in Fig. |2jb)), with suggestive, sharp transitions between low and high outputs and a 
seemingly bimodal distribution. Figs.[2jc-d) depict two more high variance segments from two 
more experimental trials. The significance of these observations may be understood from sim- 
ulated measurement signals using quantum trajectory techniques (see Appendix [B]). Fig. |2je) 
shows a simulated amplitude quadrature measurement segment, assuming perfect detection ef- 
ficiency of all photons that decay from the cavity. This simulated signal randomly switches 
between two, roughly Gaussian-distributed output states, one with a low mean and one with 
a high mean. This signal's bimodal distribution may be directly related to the bimodal, steady 
state Wigner function for the intracavity field calculated in Fig.[TJb) (discussed futher below). 
Re-simulating the same realization of the signal accounting for the calibrated overall photon 
collection efficiency (20% in the particular data set from which traces in Fig.[2ja-d) were pro- 
duced), yields a signal in which binary switching is much less distinct in Fig.[2jf), but resembles 
the single shot data shown in Figs.[2|b-d) in both the visibility and apparent time scales of the 
large fluctuations. 




Fig. 2. (a) Amplitude quadrature photocurrent data segment taken from the 14 — 19/is in- 
terval in Fig. [ijc)- The histogram of the data points in this segment are presented on the 
right, revealing a distribution that well-fits the expected normal distribution of photocur- 
rents when our cavity is empty, (b) High-variance photocurrent segment taken from the 
2 — 7/is interval in Fig.[TJc), when the atom should be near-maximally coupled to the cavity 
mode. Faint, but sharp transitions between high and low outputs and a seemingly bimodal 
distribution are visible in this single shot measurement, (c-d) Two more high variance seg- 
ments from two different experimental runs with the same parameters, (e) A typical Sfis 
amplitude quadrature segment simulated using quantum trajectory techniques, assuming 
perfect detection efficiency, with a clearly bimodal output, (f) The same simulated realiza- 
tion as (e), but with calibrated detection inefficiency yields a signal that resembles (b-d) in 
both visibility and time scale for the large fluctuations. 



To confirm the steady state quantum mechanical model in the high detection bandwidth 
regime more quantitatively, in Fig. [3] we compare amplitude- and phase-quadrature photocurrent 
distributions obtained from three different sets of mode-drive, atom-drive detunings and drive 
amplitudes (0, A, and E, respectively), with ~15% overall cavity photon detection efficiency 
(independently calibrated for each parameter set, with variations attributed to thermal drifts in 
the signal-local oscillator optical mode matching). Each histogram is computed from an aggre- 
gate of several few-/is-long segments of the highest variance photocurrent data, such as those 
presented in Figs.[2jb-d). While enhanced amplitude noise is always visible when the drive is 
first triggered, a minority of events exhibit fluctuations of near-maximal magnitude since our 
atom detection scheme is triggered by strongly coupled atoms (C = g'^/2KY± I) that may still 
be significantly less than near-maximally coupled; a maximum of C = 67 is achievable with our 
apparatus (see Appendices [A| and [B]) . Because of this post-selection to isolate near-maximally 
coupled atom transits, useful segments of photocurrent are rather sparse in our overall data set. 
We are thus limited to the presentation of small number statistics, but given the very distinctive 
features of the bistability-related signals (to be discussed below) and the straightforward nature 
of our selection criterion (highest atom-induced photocurrent variance), we are confident that 
our analysis procedure enables us to draw physically meaningful conclusions. 

For a near-detuned system at the onset of 'bistability' with {0,A,£'} = { — 1.1, .7,2.6} /c, 
wide/bimodal and narrow/normal distributions are apparent in the amplitude- and phase- 
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Fig. 3. (a-c) Wigner function representations of the expected steady state photocurrent dis- 
tributions, using cahbrated detection efficiencies and bandwidths, for three sets of detuning 
and drive parameters. Marginal distributions of homodyne measurements of any quadra- 
ture may be obtained by integrating these representations over the perpendicular quadrature 
(see Appendix [B|. (d-f) Histograms represent the phase- {cCtiI2) and amplitude-quadrature 
(Oq) photocurrent distributions from ensembles of the highest- variance segments. The his- 
tograms are compared with theoretically expected distributions obtained from correspond- 
ing Wigner functions (a-c), respectively. 

quadrature distributions, respectively (Fig. |3jd)). However, when the atom-drive detuning is 
increased to A = 3k* for the same drive amplitude, the low-amplitude transmitted field domi- 
nates (Fig.[3];e)) as the drive threshold for 'bistability' increases with |A|. The bimodal ampli- 
tude distribution reemerges when the drive amplitude is increased to £ = 3.7 /c in Fig.|3jf). 

The data are in agreement with quantum theoretical predictions (see Appendix [bJ, despite 
the use of a somewhat idealized model that assumes a static coupling rate g. Whereas g actu- 
ally depends upon atomic position and Zeeman sub-state, and can vary within a photocurrent 
segment because of complex atomic motion and imperfect optical polarization after many /is 
ll25l[7l, laser cooling and the cavity aperture ensures that atoms fall transversely through the 
mode in ^ 50/is, stochastic heating should induce diffusion <100nm during the first 5 /is of 
strong probing, and optical dipole forces should be minimal for near-maximally coupled atoms. 
While technical noises like fluctuations in the atomic coupling during a transit, cavity instabil- 
ity and laser noise could also give rise to signals with super-shot noise variance, such sources 
may be ruled out by a combination of being too small to be of importance, of occurring only on 
significantly longer time scales, and/or the fact that they would induce super shot-noise distribu- 
tions in both quadratures roughly equally due to our ^ — /c probe detuning (such quadrature- 
independent fluctuations appear to be minimal in Fig.[3j whereas the amplitude quadrature-only 
'quantum switching' is prominent). Although the complexity of the measurements' dependence 
on g precludes statistically rigorous parameter estimation, we find that a fixed effective value 
g = 0.8 X go in our analysis provides a good visual fit to the distributions of these short-time 
datasets (corresponding to a near-optimal if one makes the gross simplification of fitting the 
experimental histograms to theoretical marginal distributions with fixed g, assuming that the 




Fig. 4. The autocorrelation function for the same aggregated amphtude quadrature pho- 
tocurrent data presented as a histogram in figure [3jd) is displayed in black. Blue crosses 
represent the autocorrelation of photocurrents simulated by quantum trajectory methods 
[26 1 for identical parameters, as in Fig.[2jf). The stability of these quasi-bistable signals is 
enhanced relative to linearly scaled empty cavity transmission data taken after the atom is 
lost (red pluses, also presented with a 20MHz analog bandwidth and scaled to match the 
and 10/is autocorrelation of the high-variance 'experiment' signal), and the = 17ns 
cavity decay time (dashed orange line) characterizing the intracavity field relaxation of the 
empty resonator. We attribute the elevated and noisy autocorrelation of the atom-cavity 
transmission data at > 1/is timescales to dynamic fluctuations in g, as in Figjs] 

uncertainties of the histogram data points are constant and uncorrected) in all three parameter 
sets with ~ 10% uncertainty in g and no other free parameters. We believe that this approxima- 
tion of fixed g and the finite (20 MHz) bandwidth of our presented homodyne signals account 
for slight mismatches between theory and experiment in the amplitude quadrature splitting and 
phase-quadrature mean in the six data sets. 

Hence, even in the single-atom, ~10 photon regime, distinct high- and low-amplitude states 
of the output field are not washed out completely by quantum fluctuations. Remnant signa- 
tures of optical bistability are visible in the bimodal output photocurrents of Figs. [2] and [3] 
Similarly, as suggested in Fig.|2] we see in Fig.|4]that although the output field switches spon- 
taneously when an atom is present in the cavity, it remains correlated over timescales much 
longer than that of light transmitted through an empty cavity. This atom-induced memory ef- 
fect can be seen as an additional remnant of MBE-type optical bistability, where classically the 
high- and low-amplitude states are truly stable (due to negligible fluctuations in the system) 
and would therefore exhibit infinite correlation time. Consequently, it should also be possible 
to observe the hysteretic amplitude response characteristic of classical optical bistability by 
modulating the system drive slowly compared to the timescale for relaxation of the intracavity 
field (set in our case by the cavity decay time) but rapidly compared to the 'metastable' mem- 
ory timescale indicated in Fig. |4] Accordingly, the data in Fig. [5] were obtained by recording 
amplitude-quadrature homodyne photocurrents while sweeping the drive strength sinusoidally 
at 0.25MHz or IMHz. Figs.|5ja) and[5jb) depict representative single-shot photocurrent seg- 
ments with {0,A} = { — 1.1,.7}k' encompassing several cycles of sinusoidal drive amplitude 
modulation (AM) spanning the steady state bimodal region. Increases in both the mean and 
variance of the output photocurrent, largely in phase with the drive amplitude, can be dis- 
cerned in both of these real-time plots. However, plotting the photocurrent as a function of the 
instantaneous drive amplitude (Figs.jSjc) and[5jd)) reveals a significant hysteresis in the sys- 




Fig. 5. (a-b) Single-shot, amplitude quadrature measurements as the drive amplitude is 
swept at .25MHz and IMHz, respectively. Black traces are 20MHz bandwidth photocur- 
rents while the green, dashed traces represent the instantaneous drive amplitude, (c-d) 
Dashed red (blue) traces portray the same photocurrent data in (a) and (b), respectively, 
as a function of the instantaneous, increasing (decreasing) drive amplitude. Error bars 
represent sample mean and sample standard deviation of the same photocurrents within 
non-overlapping drive amplitude intervals. Red and blue regions represent theoretically ex- 
pected photocurrent mean and sample variance as a function of instantaneous drive (see 
Appendix [B]). (e-f) Only the sample means of 3 additional single shot measurements of 
similar duration are plotted in each figure as a function of the instantaneous drive ampli- 
tude, as in (c) and (d), overlaying the same theoretically expected photocurrent statistics. 

tern response at IMHz AM that is barely noticeable at the more adiabatic .25MHz AM rate. 
Whereas the response of the empty cavity is linear and non-hysteretic with fixed (shot-noise) 
output photocurrent variance at these modulation frequencies, nonlinear increases in the signal 
mean and variance are evident in both traces at mid-sweep. At IMHz AM, a hysteresis loop ap- 
pears to open between the upward and downward drive amplitude sweeps, with the low (high) 
state persisting over a wider range of increasing (decreasing) drive amplitudes than at .25MHz 
AM. Moreover, we emphasize that despite the large variance in the signals when viewed at 
such high bandwidth, the transition from a non-hysteretic to hysteretic mean response as the 
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Fig. 6. The essential optical and electro-optical components that define and stabilize the 
science cavity resonance frequency, drive the cQED system, detect individual atom transits, 
and perform homodyne detection of the cQED transmitted field, as explained in the text. 

AM rate increases is significant, as highlighted in Figs.[5je-f). These data are consistent with 
theoretical predictions (see AppendixjB]), again assuming the same effective value of g = .8 x go 
determined from the short- time segments analyzed in Fig.|3] 

In conclusion, we have probed the dynamic optical response of a driven, strongly coupled, 
single-atom cQED system in the vicinity of atomic saturation. While our observations further 
confirm the predictions of a fully-quantum model, qualitative remnants of semiclassical absorp- 
tive bistability are clearly visible in the data. While the quantum fluctuation-limited lifetimes of 
the low- and high- amplitude states in attojoule optical 'bistability' may be too short for direct 
use in photonic switching, our results clearly illustrate that existing theoretical models can be 
used to predict and analyze the dynamic response of real devices in the few photon regime. Such 
models can be used to make detailed predictions of the impact of quantum effects on ultra-low 
energy switch performance, which may be of interest to the nanophotonic engineering commu- 
nity. Beyond mere simulation and analysis, existing theoretical methods can be used to explore 
new approaches to the suppression of quantum fluctuations in the design of switches and related 
nanophotonic devices, for example by exploiting embedded coherent feedback control |[T0ll27ll . 
Even in a purely classical information processing paradigm, high spatial-density and ultra-low 
power nanophotonic circuit design presents intriguing new challenges for the nascent applied 
physics discipline of quantum engineering. 



Appendix 

A. Experimental apparatus 

The experiment consists of a standard cQED setup involving laser cooled ^^^Cs atoms and 
a high finesse Fabry-Perot optical resonator (J^ITll^. Through frequency and polarization 
selectivity, we attempt to drive only the {6Sii2^F = 4^mF = +4) ^ (6P3/2,F = 5,m^ = +5) 
atomic cycling transition at 852nm in order to approximate the atom as a two-level system. 
Although the cavity used in the experiment was constructed to optimize signatures of sponta- 
neous dressed-state polarization at high drive amplitudes |7 |, it proved well-suited to study the 
amplitude 'bistabiliy' dynamics presented in the main article. 

Inside a UHV 10~^ Torr) chamber and placed on a multi-stage vibration-isolation stack, 
the Fabry-Perot optical resonator is formed by two high-reflectivity (8ppm transmission, 2ppm 
loss), 10cm radius of curvature dielectric mirrors with roughly 27 /im of separation, yielding a 
300,000-finesse optical resonator for the standing wave, TEMqo, 18/im- waist transverse spatial 
mode with a field decay rate of k* = 27r x 9.3MHz. We took particular care to mount the mir- 
rors in a rotationally- symmetric manner to minimize stress-induced birefringence in the mirror 
coatings, allowing for full polarization-selectivity of the atomic transitions. The cavity length 
is tuned and actively stabilized by two shear-mode piezoelectric plates underlying the two mir- 
ror mounts. The precise cavity length and resonance frequency is continually stabilized by the 
Pound-Drever-Hall (PDH) 1291 method using an additional laser probe detuned by the desired 
probe/cavity resonance frequency by two cavity free spectral ranges (at an optical wavelength 
of roughly 826nm, which interacts negligibly with Cs). 

A Doppler-limited, magneto-optically trapped ensemble (MOT) of ~ 10^ atoms is formed 
roughly 1cm above the cavity mode in the UHV chamber. After cooling, the ensemble trap is 
switched off, allowing the cold atoms to fall under gravity towards the cavity mode and by the 
time they reach the cavity mode their free-fall velocity tends to dominate any residual thermal 
motion. Due to the strong coupling between the targeted atomic transition and the cavity mode 
(with calculated maximum value go = 27rx 56.8MHz at the cavity anti-nodes, using the dipole 
strength of the atomic transition and cavity mode volume), individual atom transits are detected 
by monitoring the (g-dependent) cavity transmission amplitude using a relatively weak and 
near-resonant probe |28|, a free space balanced photodetector, and an actively phase-locked 
optical local oscillator (LO). Although multiple atom transits per drop may be visible (apparent 
whenever C = g'^/2KY± ^1), the atomic ensemble is sufficiently diffuse such that no more 
than one strongly-coupled atom is simultaneously present in the cavity mode and we acquire 
data from only one transit per ensemble drop. It is in principle possible that a weakly coupled 
"haze" of background atoms are also coupled to the mode, but there is now a considerable 
literature (starting with 1281 ) that demonstrates that dropped-atom cavity QED systems quan- 
titatively reproduce the predictions of single-atom theory. Once a strongly coupled atom has 
been detected, the probe power and frequency shift to the desired experimental levels and data 
acquisition is initiated. 

Fig. [6] depicts many aspects of the resonance lock, optical drive, transmission measurement, 
and atom-triggering in the experiment (the laser cooling optical system is standard and remains 
tacit). A diode laser is PDH frequency locked 253MHz to the blue of a large, mechanically 
stable 'transfer' cavity mode, which is itself locked 253MHz to the red of the (65'i/25^ = 4) ^ 
(6/3/2,^ = 5), 852nm hyperfine transition (so that the diode laser is on atomic resonance). 
The purpose of this ~ 80,000 finesse, lOkHz-linewidth transfer cavity is two-fold: provide a 
stable, Cs-locked frequency reference for all the lasers and act as a cleaning cavity for the 
diode laser, producing a more narrow linewidth 852nm laser source in transmission than is 
easily achievable in the laser lock. The 852nm beam transmitted by the transfer cavity seeds 



a high power slave laser, which sources the light used in both the science cavity drive and 
optical local oscillator. The 'drive' arm of the slave output is frequency shifted to the red of 
the atomic transition again before entering an electro-optic modulator (EOM 1), which adds 
the optical carrier sideband near atomic resonance that serves as the science cavity's dynamic 
drive (the carrier and other sidebands are sufficiently detuned from the cQED system and have 
no measurable effect). The science cavity is actively stabilized by a second, titanium- sapphire 
(Ti:S) laser at roughly 826nm and frequency locked relative to the 852nm diode laser and 
the atomic transition via a frequency lock of the Ti:S laser to the transfer cavity: an optical 
sideband exactly two free spectral ranges to the red of the desired science cavity resonance near 
the atomic resonance is defined by a second EOM (EOM 2) and used to PDH resonance lock 
the science cavity in transmission. The power of the locking laser was kept at <100nW so that 
the AC Stark shift it induced in the atoms could be neglected in our analysis. As implied by 
the main article, the effects of the uncertainties in system parameters (e.g. laser noise, cavity 
resonance frequency instabilities, AC Stark shifts of the atom) are dominated by the effects of 
the variations in the atom-cavity coupling rate. 

For detecting and triggering off of strongly coupled atoms during their ~ 50/is transit through 
the mode, the drive is initially tuned several MHz from the desired frequency for the experi- 
ment and its amplitude is set below the level of atomic saturation, while the several mW optical 
LO is tuned to the experimental drive frequency. Both heterodyne quadratures of the transmit- 
ted atom-probe are detected via their interference with the optical LO. The heterodyne phase 
quadrature is used to stabilize the relative phase of the free- space LO against slow drifts in 
the signal and LO path lengths via a IkHz bandwidth optical phase lock loop. The amplitude 
quadrature is monitored by a Schmitt trigger that fires when the measured field amplitude drops 
below a threshold indicating a strongly coupled atom in the cavity. This trigger shifts the drive 
amplitude, frequency and LO-relative phase to the desired experimental configuration and ini- 
tiates data acquisition of the homodyne photocurrent at 200MS/s. Although collected at full 
bandwidth, the data presented in the main article has been additionally filtered for clarity using 
a 20MHz two-pole low pass, filtering out as much of the high frequency shotnoise as possible 
while still preserving the visibility of the fast, quantum fluctuations. 

B. cQED modeling and photocurrent predictions 

The driven Jaynes-Cummings Hamiltonian (301 ED is the standard quantum model of the in- 
ternal dynamics of a two-level atom coupled to an optical resonator mode in a frame rotating 
with the external drive (h= 1) 

H = Ag^ a -^&a^ a-^ ig{a^ a - aa^) iE{a^ - a), (1) 

where a is the cavity mode annihilation operator, a is the atomic lowering operator and ^ de- 
notes the Hermitian conjugate. From left to right, the RHS terms correspond to the atom-drive 
detuning, the mode-drive detuning, the atom-mode coupling, and the external drive, respec- 
tively. The complete quantum model of this system comes from extending these Hamiltonian 
dynamics to include processes associated with the dissipation of photons through the cavity 
mirrors at mean rate 2k per intra-cavity photon and excited atomic state spontaneous emission 
at mean rate 27^, as modeled by a and a operator-coupling to external quantum fields |[3Ql[3l1i . 
Only one of these fields, the transmitted mode, is monitored by our homodyne detection setup. 
For many applications, the entire model may be effectively represented by a master equation 
that describes the unconditional evolution of any mode and/or atom operator O (in the Heisen- 



berg picture) (321 [STl 



j^O = i[H,0\^2K(^^Oa-]^a^aO- 



1 



1 



1 



(2) 



For example, the (one dimensional) null space of the analogous, Schrodinger picture 'Liouvil- 
lian,' ^s-> corresponds to the steady state density matrix utilized in Figs.[TJb) and[3]in the main 
article. This model may be trivially extended to describe the interaction of any number of two 
level atoms with the mode (or even with multiple modes). 

Starting from this quantum model, the Max well-B loch Equations (MBEs) may be derived 
by first assuming the approximation that atom-mode operator expectations factor [13^ ^|, e.g. 
{aa'^) ^ {a) (cr^) (see also 1231 ). This approximation is equivalent to modeling the cavity mode 
as if it were a noiseless classical field coupling to a two-level atom or to an ensemble of atoms 
(as the case may be). The MBEs can take the form of a set of five real, first-order, non-linear 
differential equations of motions for {a) and other expectations. Steady state solutions for {a) 
(proportional to the mean field transmitted by a Fabry-Perot cavity) and their dynamical sta- 
bility may be found. Properly scaled, these mean-field solutions depend on the atom(s)-mode 
coupling only through a dimensionless parameter known as the 'cooperativity,' C = Ng'^/2KY±, 
where N is the number of coupled atoms and g the rate of coupling to a single atom | 6 1. Thus, as 
long as the cooperativities are equal in both cases, this model predicts equivalent steady states 
for a coupled atom and for a coupled macroscopic atomic ensemble. For highly non-linear, 
C I systems, however, the correlations between discrete excitations in the mode (photons) 
and the atom(s) that are ignored in the mean-field model may play a significant role in the over- 
all dynamics when A/^ ~ 1, as demonstrated in the main article, with our experimental C ^67, 
A^=l system. 

Measurements of the external field transmitted by the cQED system are modeled using a 
quantum stochastic differential equation (QSDE) framework |[32| |33l . Loosely speaking, this 
model represents observables of the transmitted field measured by our homodyne detection 
setup as a linear combination of cavity mode observables and 'quantum white noises.' For 
example, expected photocurrent distributions and correlations may be calculated from the char- 
acteristic functional |[33l [32l 



= (^exp psdBl,^, - P:dBo 
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where is some complex- valued scalar function of time, * denotes complex conjugation, Bout.s 
is the QSDE annihilation process of the measured, transmitted field, ^ is the time-ordering op- 
erator, and the expectation is taken over both the system and external field degrees of freedom. 

Assuming a boxcar averaging photocurrent filter of width T (so that = j3 for < ^ < T 
and Ps = otherwise), the Fourier transform of Eq. ^ produces a representation of the field 
analogous to a Wigner quasi-probability distribution 



Wt{a) 



(4) 



where pt is the reduced density matrix for the cQED system alone at time t (may be the steady 
state density matrix, for example), = P^/2Krjpta^ — P^'y/Ticrjapt is a 'measurement su- 
peroperator,' and r] is the efficiency with which photons decaying from the cavity are measured 
by the homodyne detector. As with a Wigner function, Pt{ao), the probability that a homodyne 
measurement of quadrature 6 at time t will take the (real) value ag , is obtained by integrating 



Wt{a) over the complementary quadrature, 



Pt{ae) = J dae^j,/2Wt{a). (5) 

In practice, we approximate Wt{a) by assuming short photocurrent integration times so that 
we may 'freeze' the relevant internal dynamics over the sample interval, taking 0. This 

approximation also implies the assumption that the Gaussian fluctuations in the input vacuum 
field are not correlated with system state over the effective measurement time interval. As a 
result, in this approximation Wt{a) is equivalent to a convolution of the Wigner function for 
the system state and a symmetric, mean-zero Gaussian distribution with a variance set by the 
effective integration time (which represents the contribution of vacuum field fluctuations to the 
measurement distribution). These methods and approximations were employed in Fig. [3] of the 
main article. 

Similarly, integrating Eq. ^ over some Pq^j^^2-^'^^^^^^^^^ produces the characteristic func- 
tional 

^eAk] = (^exp|/^^Ml^J,,|), (6) 

where kg is an aribtrary real- valued scalar function of time and dYg s/dt = Is is the homodyne 
photocurrent operator of quadrature 6. This functional may be used to calculate moments of 
instantaneous photocurrent measurements 1.33 >. 32 J 

(/..../.) = (-/r^^*e,.Wk.o. (7) 
For instance, the mean amplitude-quadrature photocurrent at time t is calculated to be 

{It)=Tr[-i^oPt] (8) 

where —i^opt = y/Ticrj (^pta'^ -\- apt) , while the two-time photocurrent correlation is {t <t' < 
T) 

(Itdt) = S{t -t^) -Tr[^o^xp{^{t' -t)}^oPt] . (9) 

Note that the first term on the RHS of Eq. ^ may be identified as the shotnoise contribution 
to the photocurrent correlation function, while the second is the contribution from the sys- 
tem. As the detector-filtered photocurrent operator is J^t = Jo ft {s)Isds where ft {s) is some 
filter function imposed by the detection at time t (e.g. a low-pass filter initiated at time t), the 
above two equations may in principle be used to calculate the mean and variance of filtered 
photocurrent measurements. In practice, though, we again invoke a small integration time ap- 
proximation, justified by our high bandwidth detection. If we approximate the state pt as static 
over the effective integration time (as far as detection is concerned), we may take pt in the RHS 
of Eqs. ([S]) and ^ as independent of time over the effective integration interval and approxi- 
mate ^{t^ — t) ^0, greatly simplifying these calculations. These approximations again allow 
us to model the detector filter function by a simple time-averaging filter of width T = 2/(;r/c), 
where fc = 20MHz is the cut-off frequency of the two-pole low pass filter that sets the band- 
width of the data presented in the main article. Thus the mean and variance of the photocurrent 
measurements expected from our detector when the system is in the state pt may be calculated 
from 



i^t) = tTr[-i^oPt] 
i^h = t-t^Tr[^ipt] 



(10) 



These methods and approximations were used to calculate the theory curves in Fig. [5] of the 
main article. Integrating the Schrodinger-picture master equation for AM drive (initiating sys- 
tem and drive in the ground state) returns an expected ensemble system state pt at each point in 
a drive cycle. These states were then used to calculate the expected ensemble mean and stan- 
dard deviation of the photocurrent measurements using Eqs. ([TO]), given the bandwidth of the 
data. Relating the master equation's quantum ensemble, single cycle predictions to the several 
cycle, single shot data in Fig. [5] is justified by an assumption of the system's ergodicity (see 
below). 

Finally, we relied on quantum trajectory methods |31 1, which compliment the above master 
equation-based approaches to simulate typical measurements expected from our apparatus in 
Figs. |2] and |4] While the master equation may be used to model ensembles of experimental 
realizations, a simulated quantum trajectory may be used to construct potential experimental 
homodyne measurement sequences, correctly sampled from the space of all possible sequences. 
For example, simulation of an amplitude quadrature photocurrent given a set of experimental 
parameters first involves the calculation of a possible trajectory for the internal quantum state 
vector I V^c(O) t>y numerically integrating the stochastic Schrodinger equation [26, 31] 

d\Wc{t)) -{iH^Ka^a^Y^a'^a)\\\fc{t))dt^ (^V7jc{a^a^)cdt ^dW^^^^ V2Ka\\\fc{t)) + 

[^^^{a^G^)cdt^dwP) ^^a\Wc{t)) (11) 

where {dWf^\dW^^^} are randomly generated, independent Wiener increments, {■)c denotes 
expectation with respect to |v^c(0) the state vector is forcibly re-normalized after each 
recursive update. The simulated photo-increment dYt may then be obtained using this state 
trajectory and calibrated detection efficiency r] by 

dYt = ^/fj (y2K{a + a"^ ) c + J ^ ^ ^ + 1 - r] J w/^^ (12) 

(3) 

where dWf ^ is a third, independent Wiener increment. 

Again, quantum trajectory simulations more directly reproduce time- series, experimental 
measurements, while master equation-based simulations represent ensembles of trajectories or 
measurements. In Fig. [7] we present many of the same figures in the main article using simulated 
photocurrent data based on quantum trajectories in comparison with the master equation pre- 
dictions. Figs.[7ja-b) depict simulated amplitude quadrature photocurrents with the same exper- 
imental parameters as in Figs.[TJc),[2]and[3jd) in the main article. In Fig.jTJa), the measurement 
of the the mode leaking from the cavity is simulated as detected with perfect efficiency, while 
Fig.jTJb) simulates the same realization, but with our experimentally typical 20% efficiency, as 
in Fig.[2je-f). Fig.jTJc) histograms a simulated inefficient photocurrent trajectory in comparison 
with the expected amplitude quadrature distribution using Eq. ^ (i.e. reproducing the ampli- 
tude quadrature component of Fig.[3jd) in the main article using simulated data). Finally, Fig. 
[TJd) compares the sample mean and sample standard deviation of simulated amplitude quadra- 
ture measurements during 50 cycles of a IMHz AM external drive, for comparison to Fig.jSjd) 
in the main article. Both Figs. [TJc) and (d) utilize several times more aggregate data than the 
experimental data presented in the main article in order to articulate the convergence of pre- 
dictions from quantum trajectories to those of master equation simulations. As we expect both 
types of systems to be stationary ergodic processes, we have empirically confirmed that the 
main discrepancies between quantum trajectory and steady state master equation based simula- 
tions observed in Fig. [7] arise from the marginal appropriateness of the approximation 
when modeling 20MHz bandwidth signals derived from our system. 



Simulated steady state amplitude quadrature 
Simulated steady state photocurrent with ti=1 photocurrent histogram 




Fig. 7. Reproduction of several figures in the main article replacing experimental data 
with simulated photocurrent data produced by quantum trajectory methods, (a), Simulated 
amplitude quadrature photocurrent data with {0,A,£} = { — 1.1, .7,2.6} /c, 20MHz ana- 
log bandwidth and perfect detection efficiency, (b), Same quantum trajectory realization 
as (a), but with the 20MHz bandwidth photocurrent simulated with rj = .2 efficiency, (c), 
Histogram of simulated inefficient photocurrents in comparison to the expected distribu- 
tion derived from master equation based techniques, (d), Sample mean and sample stan- 
dard deviation of simulated photocurrents from fifty IMHz AM cycles, overlaying master 
equation-based predictions. 
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